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Some years ago Gribov developed the so-called supercritical light quark confinement scenario. 

00 ■ Based on physical arguments he conjectured a drastic change in the analytical properties of the 

(«_^ ' quark propagator when the back-reaction of Goldstone bosons (pions) is considered. We investigate 

^^ , this scenario and provide numerical solutions for the quark propagator in the complex plane with 

Cn , and without the pion back-reaction. We find no evidence for the scenario Gribov advocated. As an 

^ ' aside we present a novel method to solve the quark Dyson- Schwinger equation in the complex plane 

Q and discuss new characteristics of dynamical chiral symmetry breaking in our truncation scheme. 
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I. INTRODUCTION 



The phenomenon of confinement is usually thought of as originating in the Yang-Mills 
i) ; sector of QCD. In the quenched theory with heavy sources confinement thus understood 



manifests itself in the behaviour of the Wilson loop at large distances; here an area law is 

associated with fiux-tube formation of colour-electric fields and a linear rising potential 

for heavy quarks. In the full theory, however, the colour-electric string between these 

\C) ' charges breaks due to the creation of light quark- ant iquark pairs. Therefore the potential 

oo ■ 



is no longer rising but levels out at large distances. Thus in a sense, made precise e.g. in 



cn ■ 111) 12], full QCD is not confining. 

r^ ! Nevertheless, free colour-charges are absent in the real world and the precise mechanism 

O I for this absence has to be determined in full QCD with realistic quark masses. This so- 

52 ■ called colour-confinement mechanism is still elusive even after three decades of intense 

r| ■ efforts. In a series of (partly unfinished) papers [3|, y, |5| Gribov developed a scenario 

• ^ ' of quark confinement arising from the supercriticality of colour charges. The basic idea, 

rS ; summarised in [6,, iZ], is the binding of quarks with positive kinetic energy within a bound 

c3 ' state of total negative energy. In order to guarantee a stable vacuum these states have to 

be filled up, therefore enforcing a vacuum with occupied quark states of positive kinetic 

energy in addition to the negative energy quark states of the conventional Dirac sea. 

Consequently, the Pauli principle prevents single quarks from propagating and there can 

be no corresponding asymptotic states of single quarks. 

According to Gribov [3], an essential ingredient in this picture is the appearance of 
Goldstone bosons due to the dynamical breaking of chiral symmetry. The Goldstone 
bosons, identified with the pseudoscalar pions, he conjectured to change the analytical 
structure of the quark propagator in such a way that the resulting quarks are confined by 
the supercritical mechanism. It is the purpose of this paper to critically investigate the 
actual infiuence of Goldstone bosons on this structure. 

To this end we employ a truncation scheme for the quark Dyson- Schwinger equation 
(DSE) and the quark-gluon vertex DSE developed in ref. (8|, which leads to a quark 
self-energy governed by non-perturbative gluon and pion exchange. The structure of the 



resulting DSE for the quark propagator is similar to the equation Gribov started with 
originally. In contrast to Gribov, we work with the DSEs as coupled integral equations 
rather than their - in principal equivalent - differential formulation that he favoured for 
analytical studies. However, in converting an integral equation into a tractable differential 
equation many approximations must be employed. Instead we work directly with the inte- 
gral equation and apply a truncation scheme that contains the same features implemented 
by Gribov, which has been used as a basis for hadron phenomenology and comparisons 
to lattice QCD results. In this truncation scheme we obtain information on the analytic 
structure of the quark propagator with and without pion back-reaction by a combination 
of several methods. As a result we find no evidence in favour of Gribov's conjecture. 

The paper is organised as follows. In section two we outline our truncation for the 
gluonic part of the quark-DSE, together with our approximation scheme for the hadronic 
part of the vertex, following the procedure of [^ . We specify our method for exploring 
the analytic structure of the quark-propagator in the complex plane, leaving details of the 
implementation to appendix \M In section three we present our numerical results, and 
finally give our conclusions. 

II. THE APPROXIMATION SCHEME FOR THE QUARK-DSE 
A. Gluon exchange part 

The full Dyson-Schwinger equation for the quark propagator is displayed diagrammat- 
ically in fig. [1] 
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FIG. 1: The Schwinger-Dyson equation for the fuUy dressed quark propagator 

With the dressed inverse quark propagator S^^{p) = i^A{p'^) + B{p'^) and its bare coun- 
terpart Sq^{p) = i^ + m the equation is given by 

S-\p) = Z,So\p)+9'CpZ,pJ-0^^^,S{q)T,{q,k)D^,{k), (1) 

with k = p — q, the Casimir Cp = {N^ — l)/{2Nc) and the renormalization factors Zip 
of the quark gluon vertex and Z2 of the quark propagator. In the course of this work we 
will consider the case Nf = 2 and Nc = 3. The vector and scalar dressing functions v4(p^) 
and B{p'^) can be recombined into the quark mass M{p'^) = BijP') / A{p'^) and the quark 
wave-function Zf{p'^) = l/y4(p^). The dressing functions depend on the fully dressed 
quark-gluon vertex Ti,{q, k) and the gluon propagator 

with the gluon dressing function Z{k'^). Up to, for our purposes minor, details in the 
far infrared the function Z{k'^) is well known from both lattice calculations and Dyson- 
Schwinger equations (for a review see e.g. |9|). 



Throughout this paper we will work in the Landau gauge as opposed to the choice of 
Feynman gauge Gribov adopted in his work. We expect that if the physical mechanism 
for quark confinement is triggered by the back-coupling of Goldstone bosons to the quarks 
this mechanism should be present in all continuously connected gauges. We prefer Landau 
gauge over Feynman gauge because the gluon dressing function is well known there (see 
below) and the tensor structure of the propagator is particularly simple^. 

The other input into eq. (^ is the fully dressed quark-gluon vertex T^{q, k). An approx- 
imation for the vertex in terms of the quark self-energies and the dressing function G{p'^) 



of the propagator of Faddeev- Popov ghosts has been developed in [111 . |12| . The ansatz 



ng] 



r,{q,p) = V,-''\q,p)W--''\q,p) , (3) 



with 



W-^'^\q,p) = G\{q-pf)~Z,, 
V:''\q.p) = r^^(g,p) 

2 p'^ — q-^ 

A(p2)-A(g2) 

p^ + q^ , . 

^ (p2 _ qiy + (M2(p2) + M^{q^)Y ' ^ ' 

where the ghost wave-function renormalisation Z^ has been shown to lead to a quark-DSE 
which has the correct ultraviolet asymptotic limit and respects multiplicative renormal- 
isability. In addition, the Abelian part J/f*^' of the construction is identical with the 



so-called Curtis- Pennington vertex F^ [13|. Its first three terms have been shown by 



Ball and Chiu [ij] to satisfy the Abelian Ward-Takahashi identity (WTI), 

iKT'^'^^{q,p) = S-\p)-S-\q), (5) 

with the quark momenta q and p and the gluon momentum k = q — p. As found in 
[12], the presence or absence of the scalar interaction term proportional to {p + q)u is of 
particular importance for the analytical structure of the quark propagator. In section IIIII 
we therefore contrast results obtained with the ansatz (jlj) also with the simpler vertex 



/ k k \ Z(k ) k k 

In other linear covariant gauges the propagator is given by Dfj,v(k) = ( <5^^ j^ ) — ^^^ + ^ '1.4" ; where f is the gauge 

parameter and C = 1 for Feynman gauge. Due to an exact Slavnov- Taylor identity the longitudinal part of this propagator 
remains undressed. Only for a nonlinear gauge condition can the propagator be rewritten as Z)^^(fc) = -S-«(fc^) (sl IlClll . 
where the choice of the running coupling as the dressing function is insprired from the Abelian theory. In the non- Abelian 
case such a choice already represents a combination of dressings for the gluon propagator and the quark-gluon vertex. 



which does not contain the scalar interaction term. Since the approximation ([6]) consists 
of only the first term of the Ball-Chiu solution of the WTI we will refer to it as 'IBC- 
vertex'. It represents a form of the rainbow-ladder approximation of the quark-DSE which 
has been successfully applied to the physics of light mesons |l5|]. 

We wish to emphasize, that the ansatz ([3]) for the quark-gluon vertex has similar prop- 
erties as a recent explicit solution of the quark-gluon vertex DSE [l&\. In particular the 
infrared singularity W"°-'"^''{q,p) ~ G{{q — pY) ~ (((? — p^Y'^'^ with k ~ 0.595 [1 



present in all tensor structures of the vertex is also an approximate property of the explicit 
solution which is proportional to ((g — p)^)^^/^^" [16|. We need to keep in mind, however, 
that the relative strength of the different tensor structures in the full vertex may not be 
represented well by the Curtis-Pennington part (jlj) of our vertex ansatz. This will play 
an important part in our discussion of the analytical properties of the quark propagator 
at the end of section IIIIAI 

In the quark-DSE the combination of the ghost dressing functions from the non-Abelian 
part of the vertex and the dressing function from the gluon propagator can be recombined 

to form the strong running coupling in a MOM-scheme, i.e. defined from the ghost-gluon 
vertex: 



2\ 9 r^2 I i:i\ ry / ■,'2\ 



a{e) = ^G\e)Z{k'). (7) 

The Dyson- Schwinger equation for the quark propagator then reads 

S~\p) = Z,S^\p) + CFZ,j^^^,S{q)Tt'^\q,k)P,^'^. (8) 



This equation, first developed in [llj, is quite similar to the integral equation Gribov de- 
rived his differential equation from. The quark-gluon interaction is basically given by the 
strong running coupling and the dressed vertex is chosen such that it satisfies the Abelian 
version of the Slavnov- Taylor identity. Note, however, that the present approximation is 
more sophisticated compared to the one of Gribov with respect to two points. First, the 
coupling under the integral is momentum dependent, whereas Gribov approximated even 
further by replacing a{k'^) ^ a(0). As a consequence we find the correct leading order 
anomalous dimensions for the quark dressing functions in the UV. Second, the Abelian 
part of the vertex nevertheless satisfies the full WTI as opposed to Gribov's version which 
satisfied only the differential form of the WTI valid for zero gluon momentum. For these 
reasons we believe that the approximation (jHj) is more accurate than the version of Gribov. 
The explicit expression for Z{k'^) used in this work has been determined in ref. 12 
by a fit to numerical solutions of the coupled system of DSEs for the ghost and gluon 
propagators. It is given by 

with the gluon momentum fc^, the one-loop value 7 = (— 13A'"c + 4:Nf)/{22Nc — 4:Nf) for 
the anomalous dimension of the re-summed gluon propagator and a^ = 0.2 at the renor- 
malisation scale /x^ = 170 GeV^. We use Aq^j-, = 0.5 GeV^ similar to the scale obtained 



in ref. [12] • The infrared exponent n has been determined analytically in [17|, Il8| and is 



given by K = (93 — Vl20l) /98 ^ 0.595. The running coupling a(p^) is parameterised 
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FIG. 2: The approximated Schwinger-Dyson equation for the quark propagator with efTective one-gluon exchange 
and one-pion exchange. 



such that the numerical results for Euclidean scales are accurately reproduced |12| : 
2\ _ "^(0) 47r p^ / 1 1 \ 

Here /3o = (llA^c — 2A^/)/3, and ^^(O) is the fixed point in the infrared, calculated to be 
0^5(0) = 8.915/A^c for our choice of k. Note that such a fixed point has also been found 



for the couplings from the three-gluon and four-gluon vertices |19l . |20I . 

The expressions iQ and fITOl) represent solutions of the Yang-Mills part of QCD with 
important properties. First, note the analytic structure of the gluon dressing function iQ 
produces a cut along the entire timelike fc^-axis representing the possibility of the gluon 
to decay into ghost-antighost pairs and also into gluons. However, these particles are not 
physical and need to be confined. For the gluon dressing function this is refiected in its 



spectral properties which have been determined in [1^: the gluon has a positive spectral 
function for scales below approximately one fermi, whereas it is negative for larger scales. 
As a result the gluon appears to be a free particle in perturbation theory whereas it cannot 
propagate freely at larger scales. 

The resulting running coupling (fTOj) has an analytic structure similar to the one antic- 
ipated from analytic perturbation theory [2l|. In addition, it displays an infrared fixed 
point. Thus in contrast to the setup of Gribov, where this infrared fixed point behaviour 
had to be assumed, we are in a position to use an explicitly calculated coupling with the 
same property. Note, however, that this is only possible due to our choice for the quark- 
gluon vertices (jl]) and i^. As mentioned above the explicit solution for this vertex given 
in [165 is slightly less singular than our ansatz. Using the model of [2^, which reproduces 
this behaviour, we have checked that this difference has no qualitative impact on most of 
our results with the exception of those reported in section IIIIBt where we will comment 
further. 

B. Pion back-coupling 

As stated in the introduction, Gribov argued that the effects of the back-reaction by 
the Goldstone bosons on the quarks should be crucial to generate colour- confinement 
|3|. To this end he determined a form of the pion back- reaction that couples the pion 
directly to the quark. This can be displayed diagrammatically as done in fig. [2l A similar 
expression for the pion back- reaction has been derived in ref. |8[]. Here we develop a 
modified approximation scheme leading to a slightly different interaction. 

Consider the Dyson- Schwinger equation of the fully dressed quark-gluon vertex, given 
in the first line of fig. [3l For very small momenta, a self-consistent solution to this equa- 
tion has been given in ref. 16|. Here we are primarily interested in the mid-momentum 



behavior of the vertex and in particular in hadronic contributions. To lowest order in 






+ 




+ (...) 



FIG. 3: The full, untruncated Dyson-Schwinger equation for the quark-gluon vertex [23| is shown diagrammatically 
in the first line. The second line describes the first terms of an expansion in terms of hadronic and non-hadronic 
contributions to the quark-antiquark scattering kernel. In both equations, all internal propagators are fully 
dressed. Internal dashed lines with arrows correspond to ghost propagators, curly lines to gluons and full lines to 
quark propagators. All internal propagators are fully dressed. In the second equation, the dotted line describes 
mesons, the dashed line baryons and the double lines correspond to diquarks. 
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FIG. 4: The Schwinger-Dyson equation for the quark propagator with the quark-gluon vertex from Fig. [3] (upper 
panel) and further approximated (lower panel). 



a skeleton expansion such contributions can only occur in the diagram with the bare 
quark-gluon vertex at the external gluon line. In the second line of fig. [3] we expand the 
quark-antiquark scattering amplitude of this diagram in terms of resonance contributions 
to the kernel and one-particle irreducible Green's functions. Amongst other terms dis- 
cussed in |8| one finds one-meson exchange between the quark and anti-quark lines. Of all 
the hadronic contributions this term should be dominant, since diagrams involving heavier 
mesons and baryons are suppressed by factors of ^qcd/^^h with H G {K, p,N, ...}. 

Plugging the resulting approximation for the quark-gluon vertex into the quark-DSE 
one arrives at the diagrammatic expression shown in the upper panel of fig. HI Here the 
part denoted by the subscript 'YM' denotes contributions of a purely gluonic nature. This 
Yang-Mills part of the interaction has been specified in the previous subsection. The pion- 
part is quite complicated, since it involves not only two-loop integrals but also the full 
pion Bethe-Salpeter vertex which needs to be determined from a Bethe-Salpeter equation. 

We will further simplify this expression by noting that one of the loops involves two 
bare quark-gluon vertices and a dressed gluon propagator. The latter one is suppressed 



at large momenta and at most constant if not vanishing in the infrared (see e.g. |l8l . |24 
and references therein). In our earher work [8| we have approximated this loop by the 
Bethe-Salpeter vertex, which would be justified if the full quark-gluon vertex is almost 
bare. However this leads to an over-estimation of the back-reaction and we therefore only 
assume it to be proportional to ^275''"* here. Indeed, a good agreement with lattice QCD 
results for the quark propagator and also meson phenomenology [25] is obtained by setting 
the loop to be equal to Z2J5t\ We have checked that the qualitative conclusions drawn in 
this paper do not change when employing the truncation used in ref. |8|], which is actually 
more similar to Gribov's approach. One then arrives at the approximated quark-DSE 
displayed in the lower panel of fig. HI 

Both, gluon and pion exchange are now given by a one-loop diagram with one dressed 
and one bare vertex, respectively. Compared to the previous work of ref. |8| , fig. [2l one 
of the dressed pion-quark vertices has disappeared. As an effect, the pion back-reaction 
onto the quark is somewhat reduced. This is in line with the results of [sl], where it has 
been found that the interaction of fig. [2] leads to far too strong back-reaction effects which 
finally resulted in a dramatically small pion decay constant. Furthermore, our new ap- 
proximation removes potential problems with double counting vertex contributions, which 
are generically present in DSEs with all vertices dressed. For the purpose of the present 
paper we will use the approximation of fig. HI though we also performed calculations with 
both pion-quark vertices dressed. 

In principle, the pion in the loop couples to the quark line with its full Bethe-Salpeter 
vertex function at the dressed vertex. In general this function can be decomposed into 
four different tensor structures 

Tlip, P) = r^75 {E^ip, P) - i^F^ip, P) - i^p ■ PG^ip, P) - [^ H.iP, P)} , (H) 

where r* denotes the flavour structure of the vertex, p is the relative and P the total 
momentum of the bound state. This pion bound state is the pole contribution of the 
full pseudoscalar vertex function. In the chiral limit an exact solution for the functions 
Et^,Ft^, Gt^, Ht^ in terms of the quark self energies and regular parts of the isovector axial- 
vector vertex has been given in |26| . For the leading part E.,^ of the vertex the solution in 
the chiral limit depends on the scalar part B{p'^) of the quark propagator and the pion 
decay constant /,r and is given by E^{p, P) = B{p'^)/ f^^ such that the pion vertex in this 
approximation reads 

r;(p,P) = r^75^. (12) 

For the given truncation scheme two of us have checked explicitly that this expression is 
also a very good approximation to the full amplitude E^, for a pion with realistic mass, see 
the appendix of ref. [25| . We therefore use (IT^ for the pion vertex in the second diagram 
of fig. m The resulting complete quark-DSE then reads 

S'\p) = Z,S,\p) + C,Z, j^^^,S{q)Tt'^\q,k)P,^'^ 

with k = q — p and a factor of three in the pion interaction part due to the flavour factors. 

The pion part of the quark-DSE can now be compared with the one Gribov suggested 

in ref. |3|. There are two differences. The first one concerns the appearance of two 



dressed pion-quark vertices in the back-reaction diagram considered by Gribov, whereas 
our approximation only includes one dressed vertex for reasons discussed above. The 
second concerns the form of the pion Bethe-Salpeter amplitude. Here Gribov considered 
the form 

r^(p,P)=r^{75,5-i(p)}//., (14) 

which involves further structure in the pion amplitude besides the leading 75-part con- 
sidered in eq. (1121) . Unfortunately these additional terms are proportional to the quark 
dressing function A(p^) and therefore have the wrong asymptotics at large momenta as 
compared to actual solutions of the pion Bethe-Salpeter equation. 

Finally we need to specify the values of the pion mass m^ and decay constant f^, for 
the pion propagator in the DSEs. In the chiral limit the pion is a Goldstone boson and 
r?T.7r = MeV. We use this value together with f^^ = 90 MeV. Away from the chiral limit 
we use the physical values rrin = 138 MeV and f-,, = 93 MeV for simplicity. We explicitly 
checked that the qualitative features of all our results do not depend on variations of these 
numbers. The quantitative effects are very small. 

C. Renormalisation procedure 

Before we solve equation flT3|) we have to specify our renormalisation procedure. Upon 
multiplying flT3|) with 143:4 and ^ respectively and taking the Dirac trace, one projects 
the equation onto the self-energies -B(p^) and v4(p^) contained in the fully dressed quark 
propagator S~^{p) = i'^A{p'^) + B{p'^). Schematically one obtains 

5(p2;^2) = Z2{^^^)m + Z2{|l^)TlB{p'■,^^^), (15) 

where we have made the dependence on the renormalisation point /x^ explicit. The renor- 
malisation factor Z2 is then determined by evaluating the second equation at the renor- 
malisation point, i.e. 



l + n^(/i2;;x2 



^^(/^^) = . r;.^J.,.. . (17) 



with the renormalisation condition A(/i^;/i^) = 1. In a numerical iterative procedure 
this is always the first step at every iteration step. Furthermore, away from the chiral 
limit 771 = one can eliminate the renormalisation point independent mass parameter 
m = ?7i(/i^)Zm(/i^) by subtracting (TT5l) at the renormalisation point. This results in 

B{p^; 12^) = B{^?- /i^) + Z2{^?) (nB(p'; /i') - nB(/i'; /i')) , (18) 

with the input mass B{^'^] /i^) = M{fi'^] ^j?) v4(/i^; /i^) = M^fi"^; /i^) at the renormalisation 
point /i^. 

D. Quark propagator in the complex plane 

The behaviour of the quark propagator in the complex momentum plane and the associ- 
ated analytic structure of the propagator can be investigated in two ways. One possibility 



is to read off the analytic structure from the corresponding Schwinger function 

^s,v{t) = Jd'^J^, ^''"' ^sAp') , (19) 

where as v are the scalar and the vector parts, respectively, of the dressed quark propaga- 
tor, z.e. asip'') = B{p^)/{p^A\p'') + B\^) and ay(p2) = A{p^) / {p'' A\p^) + B\p'')) . 



This method has a long history, see p^, |27|, 128|, 129| and references therein. According 
to the Osterwalder-Schrader axioms of Euclidean field theory |30|, the function asy{t) 
has to be positive to allow for asymptotic quark states in the physical sector of the state 
space of QCD. Conversely, positivity violations in the Schwinger function show that the 
corresponding asymptotic states (if present) belong to the unphysical part of the state 
space. Thus positivity violations constitute a sufficient condition for confinement. More- 
over, by fitting (Ts,v{t) with appropriate forms one obtains information on the dominant 
{i.e. closest to the origin) non-analyticity of the quark propagator in the complex plane. 
In this work we use the form 

(Ts,v{t) = boe~'''cos{b2t + h), (20) 

which corresponds to a pair of complex conjugate poles of the propagator in the timelike 
momentum plane located at rripoie = bi ±ib2. The Schwinger function is then oscillating 
around zero with periodicity 62- If &2 = one obtains an exponentially damped Schwinger 
function corresponding to a pole on the real negative momentum axis at rripoie = &i (see 
ref. jl^ for more details). Thus by calculating the Schwinger function once with and once 
without the pion back-reaction we have a reliable tool to assess possible changes in the 
analytical structure of the quark propagator. This allows us to test Gribov's conjecture. 
These findings can be further corroborated by a direct calculation of the quark prop- 
agator in the complex plane. Technically, however, there is a caveat. Consider e.g. the 
explicit form of the DSE for the scalar self-energy B{p'^) using the IBC-vertex defined in 
eq. dSD: 

V 3 r , 1 B{q^) B{{p + qY/A) 

\2n)^J ^P + m2 gM2(g2) + 52(^2)) f^ ■ ^ ) 

A similar equation holds for A(p^). Solving this equation directly at complex momenta p"^ 
entails a complex argument k"^ = {p — g)^ of the running coupling a(fc^). This coupling, 
however, has its own analytic structure, given by eq. flTUl) . Although flTUl) represents a 
reasonably motivated and justified guess for the analytic structure of the coupling we 
would rather avoid relying upon it. Fortunately it turns out that this can be easily 
accomplished. To this end we shift the loop momentum q^ in the quark-DSE such that 
the complex argument k'^ does not appear in the gluon part of the loop but instead runs 
through the internal quark part of the loop. For the shifted equation we obtain 

CfZ2 /•^4 «(?') 3i?(P) A{p')+A{e) 



B{p') = Z,m + ^jd'q 



q2 k^A^{k^) + B\k^)) 2 



V 3 f 1 Bjk^) B{{p + kf/A) 

'(27r)4y^g2 + ^2p^2(p) + 52(p)) f^ ^ y ) 
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and a similarly shifted equation for A(j9^). The running coupling is now evaluated for 
real momenta g^ only, whereas the quark propagator is determined self-consistently for 
complex momenta p^ and k"^. 

The only caveat in this procedure is the interference with the regularisation procedure. 
Obviously, if we employ a translationary invariant regularisation scheme such as dimen- 
sional regularisation such a shift would be harmless. In our numerical procedure, however, 
we use a hard cut-off scheme. Then it turns out, that the unshifted and shifted coupled 
system of quark-DSEs for A and B precisely give the same results, if Z2 is kept fixed 
while shifting. In practise we calculate Z2(yU^) for a given renormalisation point /i^ from 
the fully converged unshifted equations and plug its value into the shifted equations as 
an input parameter. Both results then agree to numerical accuracy on the real axis. We 
then use the shifted DSEs to solve for B{p^) and A{p'^) in the complex plane. Here we 
developed a new numerical algorithm, which is described in detail in appendix \K[ 

III. NUMERICAL RESULTS 
A. Numerical results on the real axis and the Schwinger function 

We first compare the quark mass function M{p^) = B{p^)/A{p'^) and the wave function 
Zf{p^) = I/AIp"^) with and without the pion back-coupling for real momenta p"^. Our 
numerical results are shown in fig. [5l In the upper panel we compare results for the 
simpler IBC- vertex (jS]), whereas in the lower panel results for the full Curtis- Pennington 
construction (jlj) are shown. In both cases we compare the solutions in the chiral limit 
and results with a small quark mass m(lOGeV) = 3MeV, which roughly corresponds to 
an up quark with m-jjg{2 GeV) = 4 MeV. For both vertex constructions the results are 
qualitatively similar. Including the pion back-reaction into the quark-DSE reduces the 
amount of dynamical chiral symmetry breaking to some extent. This reduction is larger 
in the chiral limit. Explicit values for the quark mass function at zero momentum are 
given in table [T] and agree with this observation. 

The ultraviolet behaviour of the quark mass function is given by the analytic solution 



31 



^2, _ 27rV -(^^) 

"pMllog(PVA^CD))'"'" 



M{p )asym = —^ „,,, . „,.. 7713^- + m 



tulog 



p2 



A2 



(23) 



with the anomalous dimension 7^ = ijjv^z^jv^- "^^^ quantity m is related to the current 
quark mass m in the QCD Lagrangian, whereas {^'^) is the (renormalisation point inde- 
pendent) chiral condensate. For m 7^ the dominant part of fl23l) in the far ultraviolet 
is the second logarithmic term, whereas the 1/p^-term is important at intermediate mo- 
menta. In the chiral limit this term is the only one present. From the results of fig. [5] 
we clearly infer that the condensate term, representing dynamical chiral symmetry break- 
ing, is modified by the pion back-reaction whereas the logarithmic term, representing 
the explicit breaking due to m, is not. This is in nice agreement with our expectations. 
The explicit term should be largely independent of the details of the strong interaction, 
whereas the condensate term is not. One can determine the values of the chiral condensate 
either from fitting (!23|) to the asymptotics of the quark mass function or by calculating 

- (^^), := Z,{^^') ZUf^') N^tTo J ^^.Siq'; /i^) , (24) 



11 



in the chiral limit (the trace is over Dirac indices). We determined the condensate at 
our renormahsation point /i = 10 GeV in the MOM-scheme and converted it to the 
conventional M5-result at z/ = 2 GeV using 



(^*)2GeV 



-(^^) 



MOM 
10 GeV 



/_ln(4GeVVAy 
Un(100GeVVAi,oMJ 



(25) 



with 7„ 
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llAre-2iV 



/ 



and the scales A- 



MS 



225 MeV and A 



MOM 



500 MeV. We then 



find the values given in table [T] which support our qualitative findings from fig. [5l Note 
that the unquenching effects due to the pion back-reaction are small, i. e. of the order of 
10 MeV in the (third root of the) condensate. This agrees with previous findings both in 
the DSE framework [111 . [l5| and in lattice calculations [32] . 

The effects in the quark mass function and the quark wave function are compared 
with recent lattice results of Bowman et al. 33| in fig. [61 The current quark masses 
employed on the lattice compare to m(/i^) = 16 MeV with /x = 10 GeV in our momentum 
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FIG. 5: The mass function M(p^) = _B(p^)/A(p^) of the quark and the wave function Z/(p^) with IBC-vertex 
(upper panel) and with CP-vertex (lower panel). The exphcit mass m^ is taken at /x = 10 GeV. 
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M(0) [MeV] 


(-<**)rGev)''' [^^^V] 


location of poles mripoie [MeV] 


IBC, m=0 MeV, wo % 


269 


216 


295 (5) ±i 176 (10) 


IBC, m=3 MeV, wo % 


294 


- 


322 (5) ±i 191 (10) 


IBC, m=0 MeV, with tt 


252 


208 


279 (5) ±i 160 (10) 


IBC, m=3 MeV, with tt 


278 


- 


304 (5) ±i 180 (10) 


CP, m=0 MeV, wo tt 


322 


289 


513 (10) ±i (10) 


CP, m=3 MeV, wo % 


331 


- 


530 (10) ±i (10) 


CP, m=0 MeV, with tt 


299 


276 


478 (10) ±i (10) 


CP, m=3 MeV, with tt 


309 


- 


493 (10) ±i (10) 



TABLE I: Infrared masses A/(0), chiral condensate (~(**)2(^ev) ^'^'^ P*^^® location rripoie of the resulting 
quark propagator determined from fits to the Schwinger function a"(t) for the IBC-vertex ([6]) and CP-vertex ((4| 
choice and two different bare quark masses. 
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FIG. 6: The mass function M(p^) — B{p'^)/A{p'^) of the quark and the wave function Zflp^) with CP-vertex 
compared to the lattice results of Bowman et al. 33] . The current quark masses employed on the lattice compare 
to m{n^) = 16 MeV with /x = 10 GeV in our momentum subtraction scheme. 



subtraction scheme. In the plot our result for Zfij?) is renormalised to yU = 3 GeV by a 
finite renormalisation group transformation. M{p'^) is a renormalisation group invariant. 
One finds very good qualitative and also quantitative agreement of the effects in the quark 
mass function. In the wave function we only find small effects which quantitatively agree 
with the effects on the lattice, however with a different sign. A similar difference has 
already been observed in ref. 3J| and should be clarified in future work. Apart from this 
small deviation we therefore consider the interaction defined in section Hi B I to realistically 
reproduce the pion back-reaction effects (as opposed to the stronger one considered in 
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FIG. 7: The absolute value of the Schwinger function a{t) with IBC-vertex (left diagram) and with CP-vertex 
(right diagram). 



)^. We wish to stress that the agreement of our results with the lattice data underlines 
the eligibility of our approach to critically assess the influence of Goldstone bosons on the 
analytic structure of the propagator. 

In fig. [7| we display the results from the Fourier-transform of the scalar part of the 
quark propagator. On a logarithmic scale we show the resulting Schwinger function a{t) 
for the two different choices for the quark-gluon vertex once without and once with the 
pion back-reaction. It is evident from the results that the pion back-reaction does not 
change the analytic structure of the quark propagator. Instead, it is the form of the quark- 
gluon vertex that is crucial for the form of the Schwinger function. If only the first term 
of the Ball-Chiu solution of the Abelian WTI is used, eq. ([6]) we obtain an oscillating 
Schwinger function. In the plot this oscillation is manifest in the vertical spikes. An 
ansatz representing a pair of complex conjugate poles in the momentum plane fits such 
a behaviour nicely and one obtains the pole locations reported in table [H All imaginary 
parts are clearly significant and have roughly half the size of the real parts of rUpoie- The 
inclusion of the pion back-reaction here diminishes the real parts by about 20 MeV and the 
imaginary parts by a somewhat smaller 10 MeV. The situation is totally different when 
the vertex construction (jlj) including the Curtis-Pennington ansatz is used. The relative 
strengths of the different tensor structures of the vertex are now such that the Abelian 
WTI is satisfied. From the right diagram in fig. [7] we see that the resulting Schwinger 
function now decays exponentially without any visible oscillations. Such a behaviour is 
characteristic for a quark propagator with a singularity on the real axis which may or may 
not be accompanied by a cut. The location of this singularity is also reported in table [B 
Again the inclusion of the pion back-reaction does not change the qualitative behaviour 
of the Schwinger function but merely shifts the location of the singularity by 20-30 MeV 
to lower values. 

This is the central result of this work: the inclusion of the pion back-reaction does not 



^ From the linear plot one can again clearly infer that the pion back-reaction of the quarks is largest in the infrared 
momentum region. This is in marked contrast to the findings of previous attempts to quantify the pion corrections 
[7l l35|l . We attribute this diflFcrcncc to the wrong asymptotics of the pion wave function II14I I in Gribov's equation. 
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change the analytical structure of the quark propagator. Instead, as has been discussed 
previously in [125, i^ i^ ^^^ relative strength of the scalar and vector terms in the Yang-Mills 
part of the quark-gluon interaction which is crucial. Following ref. [12!] we can interpret 
this result. Consider for a moment the Abelian theory, i.e. QED. The non-Abelian part 
y^^abei q£ ^]^g vertex is then absent and the Abelian WTl ^ exact. Gauge invariance 
then dictates a form of the fermion-photon vertex like the Curtis- Pennington construction. 
The resulting physical electron propagator then is expected to have a singularity at the 
electron mass accompanied by a cut due to the accompanying soft photon cloud. This 
agrees well with our findings. The quark-gluon vertex in the non-Abelian theory, however, 
is necessarily modified compared to the Abelian interaction. This can be seen from its 
Slavnov- Taylor identity 

G-\k^) k, T,iq, k) = S-\p) H{q,p) - Hiq,p) S-\q), (26) 

where G{k'^) is the dressing function of the Faddeev- Popov ghosts and H{q,p) the ghost- 
quark scattering kernel. It is currently not clear, though a matter of current investigations 
id . |36| whether these modifications lead to either of the singularity structures displayed in 



fig. [71 The answer to this question remains important, since oscillations in the Schwinger 
function would be a sufficient condition for quark confinement as discussed in subsection 

mm 

B. Numerical results for additional solutions of the quark-DSE 

For sufficiently strong coupling multiple solutions of the quark-DSE exist in a domain 



V = {m : < m < rricr} of the current quark- mass [37|, |38|. The appearance of multiple 
solutions is not surprising and has strong similarities with hysteresis in ferromagnets. In 
his work Gribov emphasized however that his equations for a supercritical coupling allow 



for presumably infinite different solutions in the chiral limit m ■ T his is in contrast to 



the finding in DSE investigations using the Maris- Tandy model |37l . l38l | . Here we explore 
these multiple solutions without the inclusion of the pion back-reaction, indicating that 
the different behaviour of these multiple solutions is a result of the truncation scheme 
employed and the form of the gluon interaction. 

All of these multiple solutions are connected to the perturbative running of quark mass 
at large momenta, differing only in their infrared behaviour. They can be distinguished 
by the number of zero crossings that occur in the quark self-energy B, or equivalently 
in the mass function M = B/A. The energetically preferred or physical solution is 
strictly positive definite. Without explicit symmetry breaking, i.e. in the chiral limit, the 
symmetric solution without dynamical mass generation has to exist as well. 

In fig. [8] and [9] we present various properties of the multiple solutions. In fig. [8] we follow 
the value of the mass function M(0) starting with a positive definite solution beyond rricr, 
denoted 'positive' in fig. [8] and marked with a T. Decreasing the current quark-mass to 
zero, i.e. going to the chiral limit, the positive solution is degenerate with a negative 
solution, whose mass function has the opposite sign. This is related to a Z2 symmetry 
of the quark-DSE when simultaneously sending m — >■ —m and M{p'^) —^ —M{p^). This 
symmetry is also manifest in fig. [HI Restricting our attention to M(0) positive, we now 
introduce a negative current quark mass, giving rise to a different pattern of dynamical 
mass generation and the so-called 'negative' solutions, indicated by a 'IF. 

The location of the complex poles, as obtained through the Schwinger function is shown 
in the left diagram of fig. [9] and reveals that these positive and negative solutions are 
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smoothly connected. The key difference, however, is that these negative solutions develop 
a node in the mass function as we cross m = - that is to say they are no longer positive 
semi-definite with the scalar self-energy changing sign for small {p < 10 GeV) momenta. 
On continuing to increase the magnitude of the negative mass, we reach a critical point 
nicr whose model-dependent value here is approximately 43 MeV at /i = 5 GeV for the 
IBC vertex and the interaction described in section [TTl This critical mass merely indicates 
a bifurcation where the negative solution is degenerate with the so-called Wigner solution. 
At this point, we move again towards the chiral limit, this time following a different path 
marked by a 'III'. 

Interestingly, the Wigner solution does not connect immediately to the trivial solution 
M(0) = at m = 0. Instead, we find M(0) > for m = and therefore a second non- 
trivial solution in the chiral limit. This behaviour was not observed in previous studies^ 



of these solutions [37|, |38|. There, extensive use of the Maris- Tandy [39| interaction has 

constructed to exhibit the infrared 
16] does not show the second class of 



been made. Also the soft- divergence model of |22[ 

properties of the quark-gluon vertex determined in 

non-trivial Wigner solutions observed in the right diagram of fig. [HI Both the Maris- Tandy 

model and the soft-divergence model vanish at zero momentum with different powers of 

momentum squared. The behaviour of fig. [8] is only present when the quark-gluon vertex 

is given enough strength leading to either an infrared-fixed point or a singular behaviour 

in the effective running coupling, as is the case in our interaction and that of Gribov's. 

This second Wigner solution develops a total of two zero-crossings in the mass function 
as we cross m = 0. Increasing the mass, we follow curve TV, which again bifurcates into 
two solutions at some second critical mass m ~ 0.029MeV. Continuing the procedure 
we follow path 'V and again cross m = with the now characteristic development of 
an additional node in the mass function (see fig. [9]). At this point, we make no further 
attempt to resolve such solutions since the critical mass now oscillates around m = 
with rapidly decreasing amplitude. We can only presume that such solutions continue 
to exist, with yet more nodes developing in the mass-function as the trivial solution is 
approached. We expect the multitude of these solutions to smoothly connect the location 
of the complex poles to the trivial solution, as indicated in the left diagram of fig. [91 

Once again, we point out that we did not include pion effects here, due in part to 
the ambiguity of choosing M^^ and f.^ for varying quark mass. Since Gribov found the 
same results in his approximation with the pion back-reaction [^, it seems likely that the 
inclusion of pions does not change this picture quantitatively. 

C. Numerical results in the complex plane 

To give a better picture on the influence of the pion back-reaction on the dressing 
functions at complex momentum and also to show the efficacy of the expanding shell 
method, detailed in section IIIDI and appendix [A], we show explicit solutions to the quark 
DSE in the complex-plane in fig. [TDl For the purposes of demonstration, we employed the 
IBC vertex with a small quark mass m(lOGeV) = 3 MeV and include the contribution 
from the pion back-reaction. 



^ We did, in fact, find two-nodcd solutions with the Maris- Tandy interaction. However, these have poles on the negative 
real-axis and do not smoothly connect to the positive solution as we find in fig. [5] 
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Vertex characterising the parabola [GeV^] 
Schwinger function shell method 


IBC, m=0 MeV, wo vr 


(-0.087,0) (< -0.082,0) 


IBC, m=3 MeV, wo tt 


(-0.104,0) (< -0.100,0) 


IBC, m=0 MeV, with tt 


(-0.077,0) (< -0.072,0) 


IBC, m=3 MeV, with % 


(-0.092,0) (< -0.086,0) 



TABLE II: The parabola on which the complex-conjugate poles are located, as determined via the Schwinger- 
function and by the breakdown of the direct numerical procedure. The vertex, (— m^,0) is defined as the point 
at which the parabola crosses the real-axis. 



The advantages of employing such an expanding shell method and the associated in- 
terpolation scheme becomes apparent when we use our solutions in studies of the Bethe- 
Salpeter equations. Not only are we able to provide solutions to the quark-DSE in the 
complex plane for any numerically determined gluon propagator or quark-gluon vertex 
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- without changing or making assumptions of the analytic structure by use of fit func- 
tions - but we can do this quickly and accurately. Indeed, it is gratifying to see that our 
bound-state solutions change by less than a percent for a wide selection of coarse and fine 
grids |25|]. 

As we expand our parabolas outward into the complex plane, our domain of exploration 
approaches the location of the complex poles. Close to these poles, the function becomes 
steep and cannot be reliably represented by our interpolation scheme without adaptive 
modification of the grid points. With such tuning, it is possible to see the onset of 
conjugate poles in the complex plane by looking for bumps arising in the solutions. Since 
we do not know the precise location of the pole or its residue, we cannot extend our 
parabolas any further. Moreover, there comes a point where the poles affect our numerical 
stability and lead to a breakdown of the numerical procedure. 

This breakdown allows us to determine approximately on which parabola the complex 
poles are located, where each parabola is characterised by its vertex (— m^, 0), the point at 
which the parabola crosses the real-axis. In table [III we list the vertex of the parabolas on 
which the complex pole lies, as determined from the Schwinger function and as inferred by 
the breakdown of our numerical method. Both methods are in excellent agreement. This 
provided additional justification for the Fourier-transform method exploited in subsection 

MB 

While this procedure can in principle be applied to the CP- vertex, whose solutions con- 
tain a pole on the negative real- axis, perhaps accompanied by a branch-cut, the numerical 
solution becomes much more involved due to the derivative- like terms appearing in the 
vertex. By implementing a robust numerical procedure that deals with this numerical 
singularities correctly, we believe that our method is applicable to finding solutions in the 
complex plane for such a vertex construction. We defer a detailed calculation to future 
work. 

IV. CONCLUSION 

We studied the analytic structure of the quark propagator with and without the inclu- 
sion of pion effects, in order to compare and contrast with Gribov's conjecture of quark 
confinement due to supercriticality of the colour charge. He advocated the viewpoint that 
these pions play an important role for the confinement of quarks, as indicated in the ana- 
lytic structure of their propagator. Studying a truncation scheme that essentially includes 
all features introduced and studied by Gribov, we determined the unquenching effects in 
the quark propagator due to the back-reaction of pions onto the quarks. Our numerical 
results agree nicely with corresponding lattice calculations thus underlining the reliability 
of our truncation scheme. Investigating the analytic structure of the quark propagator by 
means of its Schwinger functions and direct solutions of the quark-DSE in the complex 
plane we found that the inclusion of pion effects had no qualitative effect. This is the 
central result of our work: Gribov's conjecture does not seem to hold. 

Instead, it is the relative strengths of the various tensor components that constitute the 
fully dressed quark-gluon vertex, in particular whether these are in agreement with those 
occurring in QED due to the Ward-identity. This finding is in agreement with a previous 
investigation of the structure of the quark propagator 1^. It relegates the question of 



quark confinement due to positivity violations in its Schwinger function to a more refined 
determination of the details of the quark-gluon vertex, see e.g. ly . 
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FIG. 10: The real and imaginary parts of the quark self-energies A and B in the complex-plane for the IBC 
vertex, including the pion back-reaction. The vertical lines show the position of the complex poles, here at 
p^ = (—0.0626 ± i0.111)GeV^. Solutions without the inclusion of pion unquenching effects are qualitatively 
similar. 



We also determined the multiple solutions of the quark-DSE as a function of the current 
quark mass. In the M (p^ = 0) — -m-plane we find a behaviour similar to hysteresis effects 
in ferromagnets connected to the phenomenon of dynamical symmetry breaking. Close 
to the origin of this plane we find the critical behaviour of the Wigner solution to be 
connected to the form of the effective running coupling associated with the quark-gluon 
interaction. Employing a sufficiently strong coupling with an infrared fixed point as also 
being used in Gribov's work, we find a multitude of solutions near the trivial point. 

Our approximation scheme for the quark-gluon interaction as a composition of a Yang- 
Mills part and a part due to the pion back-reaction onto the quark is a modification of 
the one used in j3]- This modification leads to improved values for low energy constants 
as the chiral condensate and the pion decay constant and therefore has the potential to 
describe pion cloud properties of mesons and baryons via bound state equations. This is 



further explored elsewhere [25 
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FIG. 11: The bounded parabolic region in the C-plane defined by (|A2[) . for equal momentum partitioning 77 — 1/2, 
M = 1.2GeV. Solutions to the quark-DSE are required for the whole shaded region, with Re (p^) extending as far 
as the UV cut-off. 
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APPENDIX A: SOLVING THE QUARK DSE IN THE COMPLEX PLANE 



Here we give details of our numerical method for calculating the quark DSE for complex 
momenta. In Euclidean space, our quark-DSE is conventionally written as 



Sip)- 



^2'S'o + g CpZip / 



(2^ 



■lt,S{q)T^{q,k)D^^{k) 



(Al) 



with k = p — q the momentum flowing through the gluon. We want to solve this equation 
(in the previously described approximation) for complex p"^ but using real k"^ only. This 
also appears on consideration of the Bethe-Salpeter equations. For those we find ourselves 
in need of the quark propagator evaluated at momenta 



P^ 



p + TjP , p_ = p + {1 — ri)P 



(A2) 



with 1] G (0, 1) a momentum partitioning parameter and P the total momentum of the 
meson. In Euclidean space, a bound state in the rest frame has total momentum P4 = iM 
and P = 0, and so the momenta p\ of (lA2p define parabolic curves in the complex 
plane. The necessary integrals over the angles {p ■ P) leads us to require solutions to 
the quark-DSE for all complex momenta bounded by these curves. For equal momentum 
partitioning, 7] = 1/2, the region is symmetric about the real-axis (see fig. [11] for an 
example). The vertex of the parabola is located at p"^ = (— M^/4,0) and the focus at 
p2 = (0,0). 

In contemporary BS studies, one is generally forced to work in the rainbow-ladder 
approximation whereby the full quark-gluon vertex r^(g, k) of flAl|) is replaced by its bare 
counterpart 7^. We can then solve flAip for p^ G C without iteration, requiring only as 
input the quark propagator on M"*". In practice, we need only perform calculations in the 
upper half-plane, H"*", since solutions are related by complex conjugation. The caveat, 
however, is that we require knowledge of the gluon interaction for complex momenta 
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k = p — q. In general, though the analytic structure of the gluon may be surmised, the 
analytic continuation of some phenomenological ansatz is at best ill-defined. 

To avoid evaluating the gluon for complex momenta, we modify the momentum routing 
in the quark-DSE such that its momentum is manifestly real. This amounts to introducing 
a shift in the integration variable, q —>■ k, which is of course valid in any translationally 
invariant regularisation scheme. Such schemes are in general not employed in DSE stud- 
ies, for technical reasons, so one should be mindful of any boundary terms that may arise; 
with a subtractive renormalisation scheme and careful consideration of the renormalisa- 
tion conditions these spurious terms may be rigorously eliminated. What remains is the 
equation: 



Sipr' = Z,S,'+g^CFZ,F 



d^k 



YS{q)T^{q.k)D^,{k) 



(A3) 



where again we consider p^ G IHI+ along the parabolas of (]A2|) . and now k'^ G M"*", q^ G C 
Now that the integral equation depends on the a priori unknown quark propagator in the 
complex plane, we must employ an iterative scheme to obtain solutions. 

If we choose a point p^ that lies on a parabola with vertex at (— m^, 0), then the integral 
equation only requires the quark propagator be known in the region of the complex plane 
bounded by the same parabola. The most efficient way to obtain solutions is then to ex- 



pand in a series of parabolic shells stemming from the real-axis, as shown in fig. 12(b) To 



accelerate the iteration process, the previously converged shell is extrapolated outwards, 
using the Cauchy-Riemann equations, and used as an initial guess for the next shell. 

Because we are dealing with complex numbers, it is necessary to employ some 2D 
interpolation scheme. Internally, our parabolic shells are characterised by their vertex m? 



and a parameter t^, shown as stacks in fig. 12(a), which are mapped onto the parabolas 
of fig. 12(b) It is thus straightforward to take any point p^ G C and determine its 



corresponding value in (t , m ) space. Cubic-spline interpolation is used to interpolate 
along the closest two shells mf < m?' < mfj^i in t^, whilst linear interpolation in yrn? 
is sufficient for determining the value in-between. This essentially leads to interpolation 
along a parabola, such as the dashed curve shown in fig. [T21 

The drawback of this approach, however, is that without precise information about the 
location of the poles and their residues, we are unable to explore beyond the singularities 
appearing in the quark propagator. This is discussed also in the main body of this work. 
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FIG. 12: The mapping we employ to internally represent the parabolic shells and the 2D interpolation 
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